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Abstract 

(N '• 

^ \ We show how quark-hadron duality can emerge for valence spin averaged 

structure functions, and for the non-forward distributions of Deeply Virtual 
Compton Scattering. Novel factorisations of the non-forward amplitudes are 
proposed. Some implications for large angle scattering and deviations from 
, the quark counting rules are illustrated. 
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I. INTRODUCTION 

^ : 

There has been much recent interest in two related, though distinct, processes: (i) 
Deeply Virtual Compton Scattering (DVCS), measured in non-forward Compton scatter- 
^ . ing 7*(g 2 )p — > 7P (h) Bloom-Gilman duality for the imaginary part of forward 

J> ! Compton scattering, where the electroproduction of iV*'s at lower energies and momentum 
transfers empirically averages smoothly around the scaling curve F 2 (W 2 ,Q 2 ) measured at 
^ ' large momentum transfers for both proton and neutron targets [[7|. 

In this paper we shall develop a pedagogic model of the structure functions for non- 
diffractive inelastic scattering from a simple composite system that satisfies duality. Then we 
shall investigate its implications for the non-forward distribution amplitudes of DVCS |T|,rj] in 
certain kinematic regimes, compare with existing models and abstract some general features 
in hope of gaining insight into the physical significance of measurements planned at Jefferson 
Lab and HERA. 

Our point of departure is the work of Refs. |8||| , where the partons ( "quark" ) were treated 
as spinless constituents. Reference || analysed duality in the context of a large- iV c -based 
relativistic model with an empirically inspired linear potential; it focussed on the dynamics 
required for a confined struck parton to behave as though it were free, and illustrated 
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how scaling ensured. In a complementary approach, Ref. || investigated the circumstances 
whereby the structure function for inelastic scattering F(x) - whose magnitude in leading 
twist is in proportion to the (incoherent) sum of the squares of the constituent charges - 
can in general match with the excitation of individual resonances which is driven by the 
coherently summed square of constituent charges. In the present paper we shall combine 
these to show how the (non- diffract ive, "valence") scaling function emerges for a two body 
state, whose spinless constituents have arbitrary charges. 

This F(x) is intimately related to the imaginary part of the forward Compton amplitude, 
data for which lead to the probability distributions for the partons. Analogously, data on 
Deeply Virtual Compton Scattering may be used to determine non-forward parton distribu- 
tion amplitudes F(x, £, t). Here restricting ourselves to a system of two spinless constituents, 
we study the implications of the model for DVCS, and discuss the structure and physics of 
F(x,Z,t). 

The text is arranged as follows. In Sec. II, a simple model for a two-body system with 
arbitarily charged scalar constituents is introduced; structure functions and sum rules for 
the forward and non-forward Compton scattering are derived. In Sec. Ill, we show how the 
scaling functions for the forward and non-forward Compton scattering arise in this model 
and note also relations between DVCS and the forward amplitudes which may be more 
general than the specific model. Functionally the generalised parton distribution function 
for the DVCS is discussed for both Ji and Radyushkin's frames in Sec. IV. Based on the 
functionally generalised form, a phenomenological investigation of the non-forward parton 
distribution in DVCS is made in Sec. V. Summary remarks are given in Sec. VI. 



II. A SIMPLE MODEL 

The original model of Refs. |||9| ignored spin entirely and considered the inelastic scat- 
tering of a "scalar electron" via exchange of a "scalar photon" from a composite system of 
(two) spinless constituents. The differential cross section was then written as 

da " 4E/ V, (i) 



dE'dttf 16tt 2 Ei Q 4 

where the scalar coupling g had dimensions of mass, and the factor multiplying the scalar 
structure function T corresponded to the Mott cross section. 

In this paper we shall consider the more physical situation where a spin 1/2 electron 
scatters by exchanging a vector photon. The cross section then has the standard form in 
terms of two structure functions W\ and W2, which depend on v and Q 2 

<I<J ~ <tmoh[W 2 (v, Q 2 ) + 2W 1 (v, Q 2 ) tan 2 ° ' 



-[(l-e)Wi + e(l+r)W 2 ] , (2) 



O Mott 



e(l + r) 

where r = u 2 /Q 2 , and e = pll/ptt is the ratio of longitudinal and transverse density matrix 
elements. As before we consider a composite system of two spinless constituents. While this 
is still not the real world, the model contains many of the important physical features: there 
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are analogues of resonances, the inelastic structure function exhibits scale invariance and 
"quark-hadron" duality ||. 

Ref. H demonstrated the duality for the case of a single particle in a potential (effectively 
bound to an infinitely massive electrically neutral partner) . Ref. |§ extended some of those 
ideas by considering a composite state made of two equal mass scalars, i.e., "quarks" q\ and 
q 2 with charges e± and e 2 respectively at positions f\ and f 2 . It was shown how for the 
imaginary part of the forward Compton amplitude (hence the inelastic structure function), 
the excitation of resonance states of opposite parity interferes destructively in all but leading 
twist. In this paper we demonstrate how the model leads to scaling behaviour for the 
structure functions and then extend it to the imaginary part of the non-forward Compton 
amplitude (non- forward distribution function In subsequent section some general 

features for the relation between the forward and non-forward parton distribution function 
will be abstracted and discussed. 



A. Forward Compton Scattering 

The ground state wavefunction for the composite state is ^o(r), where r 1)2 = R ± r/2 
defines the centre of mass and internal spatial degrees of freedom. A photon of momentum 
q is absorbed with an amplitude proportional to J2i e i exp(z<f • fi), which excites a "resonant" 
state with angular momentum L, described by the wavefunction ip^f). 

In the previous work which ignored spin entirely 0, the inelastic structure function in 
Eq. ([[]) was given by a sum of squares of transition form factors weighted by appropriate 
kinematic factors, 

77XCLX 1 

Hv,$ = E 7pTr|%(g)| 2 <5(^ -E -v), (3) 

where q = Pi—p}, the form factors F 0i tv represent transitions from the ground state to states 
characterised by principal quantum number N(= L + 2k, where k is the radial and L the 
orbital quantum numbers). 

Focussing on the internal coordinate f, the transition amplitude is proportional to 

(e x + e 2 )[exp(iq • r/2) + exp(— iq ■ r/2)] 
+ (ei — e2)[exp(zg ■ r/2) — exp(— iq ■ r/2)]. 

The expansion, exp(iqz/2) = J2l i L Pl (cos 9) jL^qr /2)(2L + 1), projects out the even and odd 
partial waves such that the form factor is proportional to 

F ,N[L)(q) drr 2 tl)* L (r)ip (r)j L (qr/2) 

x [(ei + e 2 )5 L=even + (e 1 - e 2 )5 L=odd ] . (4) 

The resulting structure function, summed over resonance excitations, will receive contri- 
butions from L = even (odd) in proportion to (e\ ± e 2 ) 2 . Since N = L + 2k, L = even (odd) 
will imply also that N = even = 2n or N = odd = 2n + 1, and we can expose these separate 
contributions to the structure function (we factor out the charges to make them explicit 
too), 
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H^Q)= E j^[K2n(q)(ei + e 2 ) 2 + Fl 2n+1 (q)(ei-e 2 ) 2 ]5(E N -E -u) , (5) 

N(n) Ah 0^ N 



or equivalently: 



0) = E TpV-K^® + F^n+MM + e 2 2 ) + 2e ie2 (F 2 2n (g) - F 2 2n+ M)} 



- E - v) . (6) 

This simple example exposes the physics rather clearly. The excitation amplitudes to 
resonance states contain both diagonal (e 2 + e 2 ) and higher twist terms (±2eje2) in the 
flavour basis. The former set added constructively for any L and the sum over the complete 
set of states can now logically give the deep inelastic curve 0; the latter enter with opposite 
phases for even and odd L and destructively interfere. The critical feature that this exposes 
is that at least one complete set of resonances of each symmetry-type has been summed 
over. In general, the strength of the structure function will be proportional to e\ + e\ only 
if the excitation of even and odd L states sum to equal strengths. Ref. [|9]|| derived the 
circumstances under which this can occur. 

Consider scalar quarks confined in a linear potential, and described by the Klein- 
Gordon equation such that V 2 (r) = (3 A r 2 , where /3 4 is a generalized, relativistic string 
constant (compare also Ref. ||10|| ). The energy eigenvalues are then E = ±E^ where 



En = y 2(3 2 (N + 3/2) + m 2 and m is the mass of the interacting constituent ("quark"). 
This choice gives a spectrum that is in accord with that observed empirically |iT)|^T2"| since 
the energy eigenvalues follow from the similarity to the Schrodinger equation for a non- 
relativistic harmonic oscillator potential; the wave functions are algebraically as for the 
non-relativistic case and thereby enable analytic solutions. 

The contribution to J-"(q) from the N(= L + 2k) set of degenerate levels is 

from which one can immediately see that SaLo n(o) = 1- ^ * s interesting to note that 
any individual contribution, Fo } N(q), reaches its maximum value when q 2 = 2/3 2 N, at which 
point Fq N = Fq N+1 . This coincidence is true for all juxtaposed partial waves at their peaks, 
which gives a rapid approach to the equality of J2%Lo Fo2n(o) anc ^ E^=o ^0 2n+i(<z)- 

We turn now to the more physical case of real spinning electron and photon. For scat- 
tering from a system of spinless constituents, the leading contributor at large Q 2 [|13|] is the 
longitudinal response function 

R L => I J^[uW 3 -2MxW 1 ]. (8) 

The resonance sum for is analogous to that in Eq. (|^) for JF but with an extra factor of 
(E + E N ) 2 . Hence |4 



max "i 

Rl>M= E TTT-^l/0,7v(g1| 2 [(^0 + ^) 2 ^ + ^0-^)] , (9) 

N = Q 4ti Q ti N 
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where 



IfoAfil 2 = (4 + 4)[Fo,2n(Q) + Fl 2n+1 (q)\ + 2e ie2 [F 2 2n (g) - F 2 )2n+1 (g)] , (10) 

one has v max < \q\ and the sum over N denotes equivalent sum over n for N = 2n and 
N = 2n + 1. Furthermore, recall that the energy eigenvalues for the Klein-Gordon equation 



are E = ±E N where E N = ^2(3 2 {N + 3/2) + m 2 . It will be helpful now to take E N > 
and to rewrite as 



max i 

1 , . |2 



x [(E + E N ) 2 5(u + E - E N ) - (E - E N ) 2 5(u + E + E N )] . (11) 
This immediately gives the sum rule 

f + OO +°° 

S(q) = / duR L (q,u) = £ |/o^v(g)r , 

J-oo iV=0 

and hence 

/-|-CXD 
duR L (u, a) = [(el + ej) + 1e x e 2 e~^l^\ , (12) 
-oo 

which generalizes as 

/+oo 
dvR L {v, q) = [(ej + e 2 ) + 2e 1 e 2 F m {2$] . (13) 
-oo 

That this is the correct generalisation rather than F^(q) as in Ref. 0, will be obvious from 
the results of the next Section [Eq. (fTEDl. It also agrees with the physical picture in Fig. [I], 
where absorption and emission by different constituents leads to a momentum mismatch in 
the internal wavefunction of q+k — ► 2q. Ref. Q has shown how such a model satisfies scaling 
and duality in the particular case when e 2 = 0. Our results enable this to be generalised to 
constituents with arbitrary charges. We shall show this in Sec. III. First we generalise the 
above to non-forward Compton scattering amplitudes as preparation for our derivation of 
the scaling form of DVCS. 



B. Non- forward Compton Scattering 

We have calculated the imaginary part of the forward Compton scattering amplitude as 
a sum over the intermediate coherent "resonance" states. This is experimentally accessible 
by measuring d Ji^ n (eA — > eX) from the target A. We now consider the generalisation 
to the imaginary part of the non-forward amplitude: j(q)A — > 7(&)A Clearly, when we 
specialise to the case k — > q; t — > 0, we must recover our results above. 

As before, focus on the internal coordinate f, and separate the transition amplitude into 
contributions from N = 2n = even and iV = 2n + 1 = odd excited states. 

The generalization is: 
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R L (u,q,k,t) = E 



1 



„l.V) 4E E N 
N 



-(E ±E n ) 2 S(u + E tE n ) 



x 



E l4%(<f)F$n0)(ei + e 2 ) 2 + {q)^% +x {k){ei ~ e 2 ) 2 ]d L m (k ■ q) 

K L=0(1) 



where d^O) = ^/ 2Z+1 ^0(0)5 and g and A; denote the unit vector in the direction of q and 

k, respectively. Whereas in the forward case we dealt immediately with J2n \ E 0N(q)\ 2 , here 
we must first sum over the various L states within a given N. Thus (see Appendix 1) one 
now has 



q ■ k 



N 



E F oN(Q)F m (k)^j:- 

N N V P 



exp( 



f + k 2 
4/3 2 



(15) 



where due to the S function in Eq. 



( E + E )Fo N (q)F m (k) = exp 

N=even N=odd 



t = — \q — k\ 2 . Hence the sum over all states gives 

Foo(\q ~ k\) , 



{q-kf, 



A(3 2 



while that over opposite phases gives 

( E E )F 0N (q)F m (k) 

N=even N=odd 



eM-^jJ^) = F 00 (\q + k\) . 



The generalisation to the non-forward case of the sum rule Eq. (|13|) is 

F 00 (\q + k\) 



S(q, k) 



dvR L (v } q,k,t) = F 00 (\q- k\) 



[ef + e 2 2 ) + 2 ei e 2 



F 00 (\q-k\) 



(16) 



If we restrict attention to t « Q 2 , in which case \q — k\ « \q + k\, we have the sum rule 



duR L {u, q, k, t) = F oo (0(e? + ej) = F 00 {t) / dvR L {v, q) . 



This is the essence of the Ji-Radyushkin sum rule [|T],||,|T5|. 

For Compton scattering at 90°, for the case q 2 = k 2 , we have 



/+oo 
dvR L {v, q) — > (d + e 2 ) e" 
-oo 



(^/2/3 2 ) 



(ei + e 2 )^ e 1(g) 



which illustrates how the counting rules of elastic scattering at large angles 16] can emerge 



from duality. An essential part of such a derivation is the assumed mass degeneracy between 
states of a given iV but of different L. Generally, there will be L-dependent mass shifts within 
a given N level in the spectrum (e.g. in a more realistic picture with spin 1/2 qq, these can be 
induced by gluon exchange and described perturbatively with a Fermi-Breit Hamiltonian). 
These will spoil the exact L-degeneracy within any fixed N and lead to oscillations about 
a smooth s~ n dependence [0 at 90°, where s is the cm. energy squared. Such oscillations 



(14) 
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are seen in the data [17j It is important to realise that such oscillations at 90° will be a 
signal primarily for mass splitting within a given N level. This is distinct from the more 
usual measure of duality, as discussed in most of this paper, which involves the energy gap 
between multiplets of different N and different parity (and would be present even in the 
absence of gluon exchange mass splittings). 

We now establish how scaling obtains in this picture. 



III. FORWARD COMPTON SCATTERING: SCALING OF F 2 (X,Q 2 ) 

It is not immediately obvious that the above model scales. Nor do the detailed algebraic 
manipulations that follow easily reveal why this property emerges. One can begin to see 
how the scale invariance arises by focussing on the shape of the curve for an excited state 
N, and in particular the position of its maximum. This maximum is at x^j ~ l/A for A 
constituents, for any excited state N as we now show. 

The scattering is on energy shell whereby vn = En — Eq. The three momentum is 
dispersed, giving a maximum magnitude for the contribution to the structure function when 
q* N = 2(3 2 N. The energy eigenvalues E N = ^2(3 2 (N + 3/2) + m 2 satisfy E 2 N - E 2 = 2(3 2 N 
and so at the maximum we have E N — Eq = <f N) whereby (f N — v 2 N = Q 2 N = (E N — E )2E . 
Hence the peak occurs when 

Qn = 2EoVn ■ 

Thus using the kinematic variable = Q 2 /2Mi>, where M is the mass of the ground state 
hadron = J2i (i.e., summed over all constituents), then for excitation of the iV-th level, 
the peak occurs at 

(AO _ Qn __ 
Xbl ~ 2Mv N ~ M ' 

which suggests that the peak occurs at a common value of xy for all N. Consequently 
we find the physically sensible result that the peaking of the structure function occurs at 
Eq/M ~ l/A where A is the number of active constituents. To obtain the explicit form 
of the distribution we replace the sum over discrete levels by an integral, and use Stirling's 
formula (for details see Sec. V of Ref. 0, where we generalize their results to the case of 
constituents with arbitrary charge). 

Focus on Eq. (TIT]) and the leading piece at high Q 2 , i.e., the (ef + e 2 ) term. Write this, 
following Ref. |§ as 

x [(E N + E ) 2 5{E N -Eq-v) 

-(E -E N ) 2 5(E N + E + v)] , (17) 

where AiV = 1. 

Some essential steps in manipulating this are to rewrite the 5-functions, e.g. 
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6{E N - (E + u)) 



then to use Stirling's formula whereby 

1 



E 

I 2 



N 



(N+ Ej ) _ (y+E y 

{ 2B 2> 2f3 2 



1 



and to take a continuum limit J2an ~ > I Then tedious algebra gives 



(v + 2E ) 2 
R L\Q> U ) = art? 7= ex P 



u 2 + 2E u^ ( Q 2 + u 



2(3 2 



In 



v 2 + 2£ v 



Q 2 + v 2 v 2 + 2£ z/ 



2/3 2 



2/3 2 



(e? + e 2 2 ) , 
(18) 



and a similar term with v — > — v. 
As Q 2 — > oo for fixed x b j 



2Mv ' 



Q 2 



M 2 i gfl 



83^Mx bj E 
Q 2 \ 



x 







Q 2 ' 



2Mx 



2Mi 



(19) 



For spinless consituents it is Rl that has been forced to scale. The transverse response 
function Brp in this model is suppressed by 0(1/Q 2 ) [14]. This is in agreement with the 
result for spinless partons in the deep inelastic regime. The implication for the structure 
function W 2 is that 



W\W 2 (x bj , Q 2 ) = F 2 (x bj ) 



x%M 2 



e ^ U bj) (el + e 2 2 ) . 



(20) 



As discussed in Section |II A| , the higher twist term vanishes in the limit of N(n) — > 
oo, where the parity even and odd partial waves are summed to the same strengths, i.e., 

Equation (BO) immediately gives the sum rule: 



f °°%i/,.. .2 .2 

oo $bj 



~doCfoj — (3-^ ~\~ & 2 



(21) 



which can be compared to Eq. (|13D , and is analogous to the Gottfried sum rule of the naive 
parton model |1J. 

The momentum sum rule for the structure function F 2 (x b j) can be also investigated, 



r+°° 1 B 2 

J F 2 (x bj )dx bj = —(E + ^-)(e? + e 2 2 ) 



(22) 



This result sheds light on the "dynamics" of the duality picture, namely, the relation between 
resonance phenomena and naive parton model. On the one hand, it shows that the results 
of the most naive parton model, which is valid at high Q 2 and where the partons can be 
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regarded as quasi-free particles with distribution 8(x — 1/A), can be only realised at (3 — > 0. 
In this case, we find 



which is the analogue of the familiar sum rule that gives the energy fraction carried by the 
charged constituents. 

The structure function F2(xbj) exhibits an intuitive connection between resonance exci- 
tations and the naive parton model. With (3 = 0.4 GeV for the potential and m = 0.33 GeV 
for the constituent mass, we plot F 2 (xbj) versus x^ in Fig. || (solid curve); this has a broad 
and flattened behavior and is dominated by the resonance excitations. Obviously, if /3 — > 0, 
we obtain the Delta function peaking at x^j — ► 1/2 for the equal mass constituents, which 
recovers the naive parton model prediction for the valence contribution to F 2 (xbj) JO]. To 
show this, we also present the calculation (dashed curve) for a weak potential with (3 = 0.1 
GeV in contrast to the solid curve, where the dashed curve peaking at x^ = E /M is clear. 
Thus we see how the physics of the naive parton model is recovered as (3 — » 0: in this limit 
the constituents are free, all excitation levels are degenerate and the structure function col- 
lapses to an effective expression, (e\ + e 2 )5(xbj — 1/A). The physical distribution is smeared 
around x^ = 1/A\ the position of the peak is related to the number of active participants 
that share momenta, and the width of the peak is driven by their momentum spread which 
is in turn related to their confinement and the energy gap for resonance excitation. 

Equation ([£3) suggests that the resonance phenomena violate the naive parton model 
result. Recalling that the groundstate energy E = v / 3/3 2_ +7n 2 ~, where 3/3 2 ~ p% + pi 
denotes the Fermi motion momentum of a constituent, the second term in Eq. fl2"2]) can be 
understood as a kinetic energy correction to the naive parton model result associated with 
the longitudinal component of the Fermi momentum. 

Certainly, this model has many features that prevent it being a serious dynamical de- 
scription of the real world. In particular, although we have generated a pedagogic picture 
of how the sum over all the resonances leads to a scaling curve which peaks in a physically 
sensible region, the quantitative details of how the scaling function behaves away from this 
region are not sensibly produced by the model. For example, the model allows a kinematic 
range — oo < x < +oo, whereas the physical region for deep inalastic scattering is limited 
to < x < 1. 

With such caveats in mind, we now investigate the (unmeasured) non-forward DVCS 
amplitude in this model, with the goal of abstracting its general features. 



The extensive discussions of non-forward Compton scattering in the literature have used 
two frames. One (which we shall refer to as the Radyushkin's kinematics) has the momenta of 
the initial state as for forward Compton and throws all of the momentum transfer A M into the 
final state. Thus 7(g) +A(P) -> 7 (A; = q-A) + A(P + A) where t = - A 2 . Here we use A to 
denote the target system which consists of two spinless constituents. In this frame, the initial 
'j A system defines the longitudinal z direction and the final state particles have in general 
transverse components. An alternative choice is the symmetric frame of Ji, where the transfer 




(23) 



IV. NONFORWARD DVCS: SCALING OF F 2 (X, £, T) 
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A is shared by initial and final states. Thus 7(g) +A(P- A/2) -> 7^ = q- A)+A(P+ A/2). 
In this frame, both the initial and final 7 A systems have non-zero transverse momenta. 
First we define the four vectors as follows 

P=[M;0;0], 

^[^-^0;-^-,M]/2, 

A = [0; -A T ; — £M] . (24) 

For Radyushkin's kinematics one has the initial target momentum p = P and hence M = M. 
Furthermore x^ = rj ~ £. For Ji's kinematics on the other hand, one has p = (P — y) and 
M 2 = M 2 + In this frame = 

Comparison with our previous calculation is most immediate if we use Radyushkin's 
frame which has the initial jA system defining the longitudinal z direction. Results are 
identical in Ji's frame, as we shall illustrate later. 

The analogy of Eq. (|17D for nonforward scattering becomes (see also Appendix 1), 



x [{E N + E ) 2 S(E N -Eo-v)- {E - E N ) 2 S(E N + E + v)\ . (25) 

First replace k — > q — A throughout, so we consider Rl(v, q, A, t). Then after tedious 
algebraic manipulations, analogous to those used in Sec. Ill, we find a scaling behaviour as 
follows: 

Cft M 2 1 -Bo 1 9' A ^2 I A 2 n 

R L {v,q, A,t) -> F L (x bj , A, t, Q 2 ) = %&y ^ M ~ J (ej + e 2 2 ) . (26) 

Note that in this model -Foo(^) = e an d so we see the explicit presence of the elastic 
form factor multiplying a skewed distribution function where, in effect, the x^j has been 

shifted relative to the forward Compton case, in a t-dependent manner: Xbj — > Xbj{l — ^3-). 

The above was all in Radyushkin's frame. If instead we had calculated in Ji's frame, all 
the steps follow analogously leading to an identical expression to the above except that the 
argument of the exponent is modified: 




Although superficially these appear to differ, when expressed in terms of observables they 
are identical. 

To see most immediately what is happening, we first ignore corrections of 0(A 2 /Q 2 ). 
Then the argument of the exponent becomes in Radyushkin's kinematics: 

^0 . c In " 

— - X bj + t,Radl* 
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while for Ji's kinematics it is: 

M {>■ O, 

However, for the imaginary part, which is all that we are considering here, there are the 
identities [fL9| : 



= X b j/2 , 



+ 

and so in either frame and in the particular kinematics t << Q 2 , we have 

F L ~ e-^-^) 2 F 00 (t) , (27) 
which leads to the amusing factorization: 

F L (x,^,t) ^ F L (^,0,0)F d (t) , (28) 

where x = x b j- This is in contrast to a commonly used ansatz F(x, £, t) — > F(x, 0, 0)F e [(t) [|J. 
We shall return to this later at Eqs. fl5^D and (P7|). This specific example reinforces the result 
of a recent study |20| which suggested that the factorization form F(x, 0, 0)F e i(t) is not a 
general result in DVCS. 

We now show how this generalises when A 2 /Q 2 corrections are included. In this case we 
need to look more carefully at 



Q-P Qo V Q J \ Q 



^-H 1+ "VJ / l 1 -" 2 



The most general form is to write these in terms of invariants: 

q ■ A = -q ■ A = (k 2 + Q 2 - A 2 )/2 . 

The symmetry of the ensuing expressions is clearest if we define K 2 = —k 2 . A particular 
example is forward Compton scattering for which K 2 = Q 2 > 0. Although in practice we 
will be interested in K 2 = 0, we keep this variable in the formula to give 

q- A _ x bj ( i _K 2 + A 2 \ 



2q-p 2 \" Q 2 J ' 
which also confirms that rj = £ for large Q 2 . We therefore define 

m _ Q 2 fm _ K 2 + A 2 

x bj — 7, i X bj — o • 



Hence for the forward case, since t = A 2 —>■ 0, xl™ = x l b n , will be recovered when K 2 = Q 2 . 
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Then Eq. fl26|) can be manipulated into the following form 



<n fin 



t,Q 2 ) 



ef + e 



Q\K 2 + 1) 



8(3^M 



x 



M 

e w 



\ M 



-xi n ) 2 - 



x bj x lj n E 

b] bj 
& . 



2 ( E J in \1 



-(4) 



(31) 



and by steps analogous to section III, the structure function can be expressed as 

M 2 



fin 



bj-^bj 



X 



2.-1- 



/£Q, 



M 2 (xi n - 
in\2 bj 



bj> g 



6.) 



4/3- 



e 2/3 



A/ 2 j-Bo _„fin\2 



e ^ 



(32) 



which is to be compared with Eq. (p 
Thus Eq. (|32"D can be written as 



for the forward case. 







[X 



x F 2 (x 



bj ) 



-^00 (A5 



(33) 

The factor F 00 (A T ) is dependent on x^, x^ and t. When x^ — x(j n (= £) — > this form 
reduces to the parton density distribution modulated by the transverse momentum transfer 
(compare Ref. pl|). However, we can restore the F e i(t) by recognising the full import of 
Eqs. (HD and ©. 

Comparison with the forward scattering Eq. fl20|) shows that the first and third expo- 
nentials in the square bracket of Eq. ( p2|) are effectively the amplitudes for finding a parton 
at x l j™ in the initial hadron and at xlj 1 in the final hadron. From Eq. (|3^) the term outside 

the square bracket, e~^ A2//4/?2 \ is seen to be the elastic form factor for the composite state 
A, F e i(t) = F 0Q (t), which arises because the initial and final states are effectively scattered 
elastically with invariant four momentum transfer squared t = —A 2 . The central term 
in the square bracket, e M ^W^j ) ; appears to be a specific property of the sum over 
intermediate states in the dual model and merits some discussion. 

As Fig. |] schematically illustrates, we explicitly included coherent intermediate "reso- 
nant" states. Although we recover the leading twist of the quasi-free parton model, nonethe- 
less confinement is present and imposing itself throughout, in particular in the intermediate 
hadron state for which a parton has entered with xfj and departed with x{™ . If we con- 
sidered the limit (3 — > 0, recovering the most naive form of a quasi- free independent parton 
model, with no memory of confinement and no non-trivial excitation spectrum in the inter- 
mediate state, this term will vanish (unless x 1 ^ = xl™), as will the entire amplitude. As 
mentioned in section III, [3 7^ leads to a physical excitation gap, which both smears the 
distributions and leads to a finite overlap when x 1 ^ 7^ x(j n . 

It is (3 7^ that enables duality via this intermediate state overlap. Specifically, the 
basic photon-parton (Q) scattering 7(g)Q — * l(k)Q has a non-zero cross section, as does 
the hadronic process 'jA —>■ A* —>■ '-/(tyA via a specific intermediate state A*. However, 
the duality is non-trivial. If confinement is "hidden" such that [3 — > 0, then although 
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j(q)Q ~^ l(k)Q exists, its embedding in the initial and final composite states will vanish 
due to the misalignment of momenta of the struck parton relative to that of the spectator (s). 
In this case the duality is still realised because all the A* states become degenerate and 
destructively interfere unless Xyj 



x 



fin 



Thus, the structure in Eq. fl3~T|) suggests that the non-forward structure function may 
generalize to a factorisation between: 

(i) the hadron-parton distribution amplitudes; 

(ii) a "quark-in, quark-out" term associated with the intermediate coherent state; 

(iii) the invariant momentum transfer t. 

We can restore the appearance of F 00 (t) by recognising that contribution (ii) describes the 
longitudinal part of the momentum transfer being "shared" in a two-step process - excitation 
and decay of the intermediate coherent state. In effect it is an indication of so-called in- 
dependence discussed in the literature |T5| , pJ|p^ , |23|| , arising from the 
confinement in the (sum over) intermediate states. Thus Eq. 



'memory" of coherent 
may be generalised to, 



F 2 (x%,xlp,t) = JF 2 (x%)x 



1 



x F 2 (x{™) x F o(t) 



f F 2 (x%) x F 2 (x£ n ) x F 00 (A r ) . 



(34) 



Equation (|34| ) has the advantage of separating out the form factors governed by the 
external composite system [F o(t)] and internal ones. As discussed in the previous paragraph, 
such a form also explains the onset of the "quark-hadron" duality in DVCS. 

It is also interesting to note that the original non-forward structure function automati- 
cally satisfies the positivity constraint of Ref. [231 , 



In Eq. ( PS] ) , one can treat the non-forward 
structure function as a product of two forward ones with incoming photon momenta q and k 
associated with a rotation function. Apart from the rotation function, the positivity bound 
for the non-forward structure function thus can be realized. Equation ( |34"D possesses the 
same feature, in which the positivity bound gives an inequality as follows: 



F 2 {x%x$,t)< y/F 2 (a^)xF 2 (x^) 



(35) 



An alternative way to write Eq. (]5ID is 



F L {xT v x f C^Q 2 ) 



(el + el) 



Q 2 (K 2 + t) 



8f3^MJx%x f b f l E 



fin 



M 2 I Bp ^bj^^bj 

~yr \ M 2 



(36) 



which can be generalised to 

F L (x*",xg n ,t,g 2 ) = F L (x,t } Q 2 ) x F 00 (t) x 



(37) 



where we have defined: 



13 



^(< + 4T)/2; £ (.r;;;-,-(;")/2. (38) 

Accordingly, in this case the structure function can be written as 

F 2 {x,t,t) = (e 2 1 + ei) {x ~ [ X + F 2 (x)F el (t) , (39) 

x z 

where F2(x) is the structure function for the forward scattering. 

The reduction of Eq. ( j39|) to forward Compton scattering is obvious. If there were no 
intermediate coherent state, we expect that the physics would be symmetric in and £&j n , 

J bj — ^bj 



through which the forward Compton scattering (xfj = £jy) corresponds to the kinematics, 
x = x 1 ^ = x(j U and £ — > 0. We conjecture that if the intermediate coherent state were 
absent, as in the incoherent parton model description, a plausible factorisation between the 
overall F el (t) and parton effective probability distribution, symmetric in x 1 ^ and x{™ would 
be 

= 4P,t) = F 2 (x,£ = 0,t) = F 2 (x) x F el (t) . (40) 

Note that x = (x|," + xjy )/2 and so this contains Eq. (p8[) as a special case in the approxi- 
mations that were used there, whereby x(j n — > and, as remarked earlier, differing from a 
commonly used ansatz. 

In the Gaussian wavefunction model, the equivalence between Eq. (|3"4D and fl3Tf) is exact. 



In Eq. ( [TOD the non-forward amplitude effectively becomes an "average" of the forward 
distributions, modulated by the kinematic factor [1 — (£/x) 2 ] where £ — > in the forward case. 
It remains to be investigated whether this analytic equivalence between Eqs. ( |34"D and (p9[) is 
an artifact of the Gaussian distributions and the restriction to Fl with spinless constituents. 
The latter in particular needs study as the imaginary part of DVCS is measurable in electron 
scattering from a polarised target, and with spinless constituents this is manifestly a non- 
leading effect. 



V. PHENOMENOLOGY 

In this paper we have constructed a model for composite systems consisting of two spinless 
particles. The functional relations between the forward and non-forward distributions, and 
the elastic form factor of the composite system, are summarised in Eqs. ( |34"D and (j39[). 
For phenomenology we want to investigate the implications of these factorisations for the 
non-diffractive structure functions of the nucleon. We adopt a dipole elastic form factor 
Fait) = 1/(1 — t/ty 2 , where A = 0.7 (GeV/c) 2 is the empirical energy scale, and choose 
\t\ = 1 (GeV/c) 2 , whereby the factorisation of Eq. (0) implies that the non-forward parton 
distribution will have the form: 

F 2 fe, I/ , ( ) = V^^ TT (^|^- («) 

Moreover, for more realistic predictions relevant to a three-body constituent system, we 
adopt the commonly used parametrisation F 2 (x) = x(l — x) 3 . The results are displayed in 
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Fig. |. As Xi(xf) < 1, one has (xi — Xf) 2 « 1 and the F 2 (xi,Xf,t) is almost symmetric 
in Xi and Xf. Hence for fixed Xf the figure shows a rise and fall as a function of Xi similar 
to that observed for the "usual" forward structure function, and to a bag model study by 



Ji et al |15|. However, there are significant dynamical differences with that model: we find 



the distribution exhibits a sensitivity to Xf, which is not apparent in In contrast to 



Ref. [15], in which the origins of the £ insensitivity are obscure | 25fl , the systematics of the 
distribution in the present model, especially its Xf dependence, are more transparent. 

For comparison, we also investigate the factorisation of Eq. (|39|) using the same quark 
counting rule parametrization as input. The original form only applies in the limit > x^ 
and so it is not clear whether its generalisation should be taken seriously out of that region. 
However, if we take the generalized form to apply throughout the entire physical range where 
the variable x and £ each have — 1 < x < 1 and £ > 0, we obtain some interesting features. As 
shown in Fig. ^, with \t\ = 1 (GeV/c) 2 , the distribution function is positive at x > £, which 
is the region dominated by the constituent (quark) distributions. The distribution function 
exhibits significant sensitivity to £. The region x >> £, which corresponds to Q 2 ~ \t\, is 
the place where resonance effects could play a role if \t\ is small. As is to be expected with 
our quark counting rule parameterization, the largest probabilities of the quark distribution 
occur at ~ 1/3. A sign change occurs at x — £, where Q 2 » \t\. In the region x < £ 
where K 2 +t > Q 2 , the partons have negative momentum fraction x{* n = x — £ < 0; this 
may be interpreted as due to antiquark q which can play a role in this kinematic region. 
Interestingly, the interpretation of this generalisation is consistent with that in Ji's frame 
(see e.g. discussions in Ref. [0]). One can see that with £ — > 0, which corresponds to larger 
Q 2 with a fixed t, this factorisation succeeds in reproducing the forward Compton scattering. 
A scaling behavior is also observed. However, notice that a fast crossover occurs at x = 
if x < £, where the usual quark density becomes infinite. This feature could imply that the 
physical region of this generalisation form only make sense at x > £ p5|,p6| . 



VI. CONCLUSION AND DISCUSSION 

We have an explicit model, in which both forward and non-forward Compton scattering 
have been investigated. The model also exhibits scaling. Although this model is far from 
reality, we may still draw intuition from it as to how the physical DVCS and related processes 
may behave. 

For forward scattering we have extended the work of Ref. M to constituents with arbi- 
trary charges and given an explicit dynamical model realising the general picture outlined by 
Close and Isgur ||. For the first time we have also determined the implications of such mod- 
els for the non-forward Comption scattering. We stress that this is at best only a pedagogic 
picture due to its emphasis on spinless constituents. 

The model showed scaling behavior explicitly for the forward scattering and also satisfied 
well known sum rules. This also has been extended to the non-forward case. We find here 
too that scaling is predicted and explicit £ dependence is also seen. 

The behaviour of the Compton amplitude at 90° showed how the effective simple s~ n 
dependence, previously derived from counting rules |Tj|, arises. For us it is directly driven 



by the dependence of F e i(t) once our Gaussian forms are generalised to phenomenological 
form factors. We found that the degeneracy among states with a common N(= L + 2k) but 
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different L causes a destructive interference among all but the elastic Born term. Thus the 
s~ n behaviour, in this interpretation, is effectively dominated by elastic scattering in the 
direct channel. Away from 90° in the model this interference is no longer exact. Indeed, in 
the real world where this N-L degeneracy is broken (e.g. by spin dependent effects arising 
from one-gluon exchange), oscillations around the overall average smooth s~ n are predicted 
rather naturally. It is therefore interesting that such effects are qualitatively evident in the 
data |T7||. A quantitative investigation of this will be reported elsewhere. 

Our work is restricted to the imaginary part of DVCS. This in turn is measurable in 
electron scattering with polarized target. This exposes the limit of our spinless model. It is 
the longitudinal response functions that are leading in such a model and exhibit interesting 
factorisations (e.g. Eq. fl34[) and fl39"l)). We expect that such relations will survive for 
investigations with more realistic models. However, it is perhaps interesting to note that 
the factorizations exhibited in the present model are novel and different from some popular 
ansatz used in the literature M. 
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APPENDIX 1 

For a given excited level N, there are degenerate states with L = N, N — 2, • • • 0(1) for 
N = even(odd). Whereas in the forward scattering we dealt immediately with J2n \Fon(q)\ 2 , 
here we must first sum over the various L states within a given N, thus 

E f E 4i\k)^mUo) 

N \L=0(1) 

where 9 is the relative angle between the incoming and outgoing photon momentum, i.e. 
k-q = kq cos(#) (in what follows q and k denote the magnitudes of the three- vector momenta q 
and k, respectively). In the potential model of Ref. |J, after sum over L the above expression 
reduces to 

w) ^'-^ (42) 

As illustration we show the first non-trivial case for N = 2. Here, for L = 2 
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and for L = 
Thus 

F 02 (q)F 20 (k) = F$(k)F$(q)d 2 00 (9) + J® ®t&{0) , 

which is the N = 2 component in Eq. (f42|). 

That this result generalises is readily seen in the Cartesian basis where the sum over 
resonances for non-forward Compton scattering in this relativistic model. We start with the 
general form 

M = 5(q + Pi-k- P f ) £<*o(r)|[eie- ife> / 2 + e 2 e l ^ 2 ]\^ N ) 

N 

x(^|[ ei e^/ 2 + e 2 e-^ 2 ]|^ (r)) , (43) 

where = VwVvV'riz is the harmonic oscillator wave function, and N = n x + n y + n z is 
the main quantum number. The one dimension harmonic oscillator, e.g. the z component 
has the expression 

1/2 

e-^ 2 H n XPz/2) , (44) 
where H nz (f3z/2) is the Hermite polynomial and the orthogonal relation is, 

/+oo , 
H n (x)H m (x)e- x dx = V^2 n n\Snm ■ (45) 
-oo 

For the excitation process we fix the incoming photon momentum to be colinear with that 
of the target, and refer to it here as the z direction. Only the z compoment ip nz will be 
excited. We have 

faM^M*)) = ~^ ( j^f e~^ 2 . (46) 

— * 

The excited intermediate state will emit a photon with momentum k and fall to the 
ground state. The transition amplitude thus can be expressed as 

(^o(x)^o(l/)^o(^)|e- j(fc ^ + ^ +fc ^ )/ Vo(x)^o(y)^(^)) • (47) 
Following the same strategy, we can explicitly derive, 

(Mz)\e- lk > z/2 \^nM) = ^ {-^Y > ( 48 ) 

for the z component. The x and y components are 







2^2 n *n z 
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(My)\e- lkyy/2 \My)) = e- k2 ^ 2 . 

Thus, we obtain 

(Mx)My)M*)\^ /2 \Mx)Mv)1>n,{z)) - 



(49) 



ik cos 6* x 



-fc 2 /4/3 2 



v/2/3 



(50) 



where fc 2 = kcosd, and # is the angle between the k and z direction (the direction of q). 
The total transition of Eq. (B3f) can be then expressed as 



M = 5(q + Pi-k-P f ) £ e 

^=0 



-(fc 2 +,f)/4/3 2 



X 



. 9 9 , 1 / ikq cos 9\ N „ 1 / ikq cos 9\ N 



(51) 



where the main quantum number iV = n 2 in this transition (which in general will include 
degenerate states with L = N, N — 2 • • •). The above deduction gives the origin of Eqs. ([42]) 
and ©. 



APPENDIX 2 



In the previous Sections, the term proportional to eie 2 has not been discussed in detail. 
Here, we shall show that this term would vanish in the limit of N —>■ oo for the non-forward 
Compton scattering process. 

The positive energy term proportional to e\e 2 in the structure function (Eq. @) is, 



N 



N- 



(52) 



Correspondingly, the term proportional to e\e 2 in the scaling function is, 



Jfe(„, g, fc, t) = 2e ie2 £ ^ 2 M -2^ 



2V 



e ~^5(E N - E -v) 



(53) 



where |g| 2 = 2/3 2 N. Compared to the term proportional to (ef + e 2 .), a factor (— 1)^ arises 
from the sum over the resonances. We therefore separately consider the iV = even and odd 
terms which gives, 



R L (u,q,k,t)=Rr n + R° L M 

<r 2 +f 

= 2e 1 e 2 e 



2/i 



^W^TW 1 q-k 



E 

n=0 



n=0 



4E f3 2 (2n)\ \2(3 2 

2n+l 

S(E 2n +i -Eq-v) 



q- k 

4E (3 2 (2n + l)! I 2/J 



(54) 
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Following the analogous steps as outlined for the term of (ef + e|) in Sec. Ill, the factorial 
part can be expanded the same way by using the Stirling's formula. After some tedious, but 
essentially the same algebra, we obtain, 

j^even _|_ j^odd 

x [5(E 2n — Eq — v) — S(E 2n+1 -E - u)\ , (55) 

where the 5 functions are shown explicitly because we assume that the photon energy v 
always satisfies the condition that the excited states are 2n for R e ^ en and 2n + 1 for R°^ d . 
Consequently, taking the limit n — » oo, E2 n = -^2n+i can be satisfied, which thus leads to 
j^even _|_ j^odd _ q Namely, the twist term vanishes. Similar investigation for the energy 
negative solution results in the same conclusion. 
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FIG. 3. Structure function F 2 (xbj) without the charge factor. (5 denotes the linear potential 
strength. 




FIG. 4. Schematic diagram for the quark-hadron duality in the non-forward Compton scat- 
tering. 
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FIG. 5. Non-forward structure function following the phenomenology of Eq. (f4l|) with t 
(GeV/c) 2 fixed. As Xf/xi ~ —t/Q 2 , Q 2 increases with increasing ccj. 




FIG. 6. Non- forward structure function following factorisation of Eq. (p9|) with t = — 1 
(GeV/c) 2 fixed. Note that t -> as x £. 
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